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Abstract 

Background: The genetic diversity of Mycoplasma agalactiae (MA) isolates collected in Spain from goats in an area 
with contagious agalactia (CA) was assessed using a set of validated and new molecular typing methods. Validated 
methods included pulsed field gel electrophoresis (PFGE), variable number of tandem repeats (VNTR) typing, and 
Southern blot hybridization using a set of MA DNA probes, including those for typing the vpma genes repertoire. 
New approaches were based on PCR and targeted genomic regions that diverged between strains as defined by in 
silico genomic comparisons of sequenced MA genomes. 

Results: Overall, the data showed that all typing tools yielded consistent results, with the VNTR analyses being the 
most rapid method to differentiate the MA isolates with a discriminatory ability comparable to that of PFGE and of 
a set of new PCR assays. All molecular typing approaches indicated that the Spanish isolates from the endemic area 
in Murcia were very diverse, with different clonal isolates probably restricted to separate, but geographically close, 
local areas. 

Conclusions: The important genetic diversity of MA observed in infected goats from Spain contrasts with the 
overall homogeneity of the genomic background encountered in MA from sheep with CA in Southern France or 
Italy, suggesting that assessment of the disease status in endemic areas may require different approaches in sheep 
and in goats. A number of congruent sub-typing tools are now available for the differentiation of caprine isolates 
with comparable discriminatory powers. 
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Background 

Contagious agalactia (CA) is a disease of small rumi- 
nants associated with several clinical signs, including 
mastitis, arthritis, keratoconjunctivitis, pneumonia and 
septicaemia. CA is listed by the World Organization for 
Animal Health and is responsible for significant eco- 
nomic loss, with its severity and expression being 
dependent on several factors, including the host species, 
the aetiological agent, the production system and the en- 
vironmental conditions [1-4]. In this context, defining 
the sanitary status of herds or regions with respect to 
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CA is essential, but also very challenging. Three statuses 
are currently recognised: disease-free areas, areas with 
sporadic acute outbreaks and endemic areas, where in- 
fection is widespread although not always acute [1,4]. In 
CA endemic regions animals typically have no to transi- 
tory clinical sign and there are many asymptomatic aur- 
icular carriers [5-8]. Preventive and therapeutic 
strategies remain very inefficient in the control of CA, 
most likely because of both pathogen-specific features 
and the lack of epidemiological data. Thus far manage- 
ment strategies have been the most satisfactory method 
for controlling this disease [4]. 

Mycoplasma species responsible for CA include Myco- 
plasma mycoides subspecies capri (Mmc), Mycoplasma 
capricolum subspecies capricolum (Mcc) and Mycoplasma 
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Table 1 Molecular typing of the Mycoplasma agalactiae isolates and strains 
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a Code for the geographical origin of the isolate and the year of isolation: YE: Yecla, MA: Malaga, LO: Lorca, ML): Murcia, JU: Jumilla. The symbol < is used when the exact year is unknown (Unk). 
b Profiles as defined by PFGE analyses. 

c Number of repeats based on the consensus sequence determined for each VNTR. * indicates VNTR for which the PCR product was sequenced. 
d Profiles as defined by VNTR analyses based on data a combination of data generated with each VNTR. 
e vpma profiles defined based on data in Additional file 2: Table S2. 
f AGA-BOV profiles defined based on data in Additional file 2: Table S2. 

9 PG2: data obtained were identical to those obtained with strain PG2; 5632: data obtained were identical to those obtained with strain 5632. 
BTM Bulk tank milk. 
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putrefaciens, from the spiroplasma phylogenetic group, 
and Mycoplasma agalactiae (MA), from the hominis 
phylogenetic group, with this latter species being desig- 
nated the stricto sensu aetiological agent of CA [1,4]. 
Interestingly, CA in sheep is often associated only with 
MA, while the situation can be more complex in goats, 
in which CA is commonly associated with both MA and 
mycoplasmas of the "mycoides" cluster [1,4,8]. 

The complex clinical situation described above 
requires tools capable of exploring the dynamics of in- 
fection in order to set up management strategies. Thus 
far, several typing techniques have been used to 
characterize CA agents, including analysis of variable 
numbers of tandem repeats (VNTR) [9], multilocus se- 
quence typing (MLST) [10], Southern blots using DNA 
probes [11], insertion sequence (IS) typing [8,12] and 
pulsed field gel electrophoresis (PFGE) [8,13,14]. How- 
ever, their capacity to answer epidemiological questions 
is variable. For instance, PFGE is very useful when com- 
paring an asymptomatic goat herd infected with multiple 
polymorphic Mmc strains to an outbreak in which one 
unique clone was predominant [8]. However, application 
of the same technique to MA isolates in another situ- 
ation resulted in poor discrimination and was not in- 
formative about the strain circulation [13]. 

A number of studies has been conducted to 
characterize and compare isolates of MA from outbreaks 
of CA. These have generally focused on sets of isolates 
from various geographical origins [11,12,14], but the 
reports have provided few details about the status of the 
herds with respect to infection, the type of herds or the 
animals from which the isolates were collected [9,10]. In 
addition, typing analyses that have been performed in 
well-documented epidemiological contexts have mainly 
examined CA in sheep [13,15] leaving the situation in 
goats unexplored. This gap is surprising, as CA in goats 
can be of considerable economical importance, and also 
because the two currently available genomic sequences 
[16,17] from which many typing tools were derived, are 
both from strains with a caprine origin. 

The present study is, to our knowledge, the first to con- 
duct a detailed characterization of potential molecular dif- 
ferences between caprine MA isolates collected in a 
geographically restricted area (Murcia, Spain). This area 
hosts a native goat breed, the Murciano Granadina, and 
uses a semi-intensive production system based around a 
single insemination centre. Although this limits the entry 
of animals from other regions, Murcia is known to have 
had endemic CA for many years [18,19]. This has raised 
the question of the mode of dissemination and the origins 
of the agent responsible for CA in this area. 

The purpose of the present study was to assess the rele- 
vance of typing approaches already validated on large and 
diverse strain collections (PFGE, VNTR, IS-typing) to 
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Figure 1 PFGE showing representative Sma\ profiles. Lanes 1 
and 2: Reference strains PG2 and 5632, respectively; lanes 3 and 4: 
isolates AG4 and AG28 (profile A); lane 5: isolate AG18 (profile A1); 
lane 6: isolate AG26 (profile B); lane 7: isolate AG14 (profile C). 
Molecular sizes (X ladder, Bio-Rad) are indicated on the left in kb. 
The lowest, detectable, PG2 band (Lanel) is predicted in silico to be 
of approximately 9 kpb. 



discriminate a limited number of caprine isolates from a 
common geographical origin (Murcia) and to explore the 
discriminative power of new, more specific markers. This 
is an essential step before addressing more complex epi- 
demiological questions that can be crucial in improving 
management of CA in goats in endemic areas. 

Results 

Strain typing by PFGE 

Four restriction endonucleases, Sma\, Mlul, Kpnl and 
BstTl, were selected that yielded easily interpretable 
PFGE banding patterns on our set of isolates. Their pro- 
files clustered the 13 field isolates in 3 main groups A, 
B, C (Table 1), as illustrated by a representative gel in 
Figure 1. Group A included 6 of the 13 isolates, as well 
as the PG2 type strain. This group was very similar to 
the main European MA profile reported in 2008 by 



De la Fe et al. BMC Veterinary Research 2012, 8:146 
http://www.biomedcentral.eom/1746-6148/8/146 



Page 4 of 9 



Table 2 Features of VNTRs Mycoplasma 
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VNTR Unit length (in nt) 
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PG2 
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CTATTATTAC 






AG32 
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CTCTATTAA 






5632 


6.4 
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TTACCTTT 




19 21 


PG2 


4.9 
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TGI I I ICI IGC 
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AG14 


4.9 


186 


TTCTTCTTGT 





Proportion of sequence conserved between the strains examined. 



McAuliffe et al. [9]. Interestingly, one isolate (AG 18) 
yielded a pattern that differed slightly from that of the 
other members of the group. After Smal digestion, it 
yielded the 6 characteristic fragments as well as a smal- 
ler one of 45 kb (Figure 1). Analysis with other endonu- 
cleases confirmed that this isolate appeared to have a 
genome that was larger than the other group A isolates 
by 40 to 50 kb. Group B included all isolates from bulk 
tank milk (BTM) samples collected in Jumilla between 
the months of January and March 2000. These isolates 
appeared to be strictly clonal. Group C included two iso- 
lates from semen that were collected in 2008 from two 
different goat males in Lorca. No field strain had the 
same restriction pattern as the sequenced 5632 strain. 

Previous PFGE analyses conducted on MA proved 
to discriminate poorly between local isolates from Italy 
[13] and, more recently, between isolates from all over 
Europe, all of each were found to be very similar to PG2 
[9] . However, in the present study PFGE pulsotypes were 
not homogeneous, despite the restricted sampling area, 
and the three groups of strains appeared to be distinct 
with a discrimination index (D) of 0.73 according to 
Hunter and Gaston [20]. 

Correlation between PFGE profiles and VNTR analyses 

Strains were subjected to VNTR analysis targeting 4 
VNTRs, designated VNTR5, 14, 17, and 19, as previ- 
ously described [9]. PCR results were first analysed by 
gel electrophoresis, which detected no difference be- 
tween strains for VNTR5 or VNTR19, with all profiles 
identical, including those of the PG2 type strain and of 
5632. While VNTR14 was of limited value as it only 
separated the strains into two groups, results obtained 
with VNTR 17 divided the 13 field isolates into 3 main 
clusters (TR2, TR3 and TR4) that matched the clusters 
obtained when comparing PFGE profiles (Table 1). 



These data were further supported by comparing 
sequences of the PCR products obtained from represen- 
tatives of the different VNTR profiles (Table 2). For this 
purpose, PCR products from (i) VNTR 5 of isolates 
AG13, AG28 and PG2, (ii) VNTR 14 of isolates AG4 
and PG2 (iii) VNTR 17 of PG2 (PI), AG30 (P2), AG32 
(P3) and 5632 (P4) (Figure 2), and (iv) VNTR 19 of 
AG 14 and PG2 were directly sequenced (see Additional 
file 1: Table SI for each VNTR sequence). Based on 
these sequence data, a consensus sequence was gener- 
ated (Figure 2C for VNTR17) and the number of tandem 
repeats defined for each VNTR (Table 2). A VNTR pro- 
file, TR1 to TR5 (Table 1), was ascribed to each strain 
and the overall VNTR analysis found to match the PFGE 
clustering, although it did not discriminate strain AG 18, 
which had a pulsotype slightly different from the overall 
group A. This is reflected by the index of discrimination 
which is slightly lower for VNTR typing (0.67) than for 
PFGE (0.73). 

Using the vpma profiles as a typing tool 

In MA, the vpma locus has been thoroughly studied and 
shown to encode genes with related sequences [15,21- 
23]. The size of the vpma repertoire (number of vpma 
genes) varies between strains, but is constant within 
each strain. Here, we assessed the vpma profiles of our 
set of strains by Southern blot analyses using genomic 
DNA digested with ^4sel or Alwl and hybridised to dif- 
ferent probes. The first probe was a single digoxigenin- 
labelled oligonucleotide probe, A3F, that detects the 
conserved 5' ends of all vpma genes and thus assesses 
(when appropriate restriction endonucleases are used) 
the number of vpma genes in the strain, regardless of 
their specific vpma sequences. In contrast, the other 6 
probes individually target each vpma gene found in the 
PG2 type strain {vpmaV to vpmaZ). Based on the 
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Figure 2 Profiles obtained with (A) VNTR 17, (B) VNTR 5 and (C) VNTR 17 nucleotide repeat sequences corresponding to PCR products 
of individual profile (PO to P4). Positions with similar nucleotides in all profiles are marked with an asterisk. Isolates are indicated above the first 
gel. Lane 1, size marker (Smart Ladder Eurogentec). 



hybridization results, a vpma profile was assigned to each 
strain (Table 1 and Additional file 2: Table S2 for detailed 
results). 

As shown in Table 1, the vpma profiles were slightly 
more discriminative (D = 0.74) than the VNTR or PFGE 
profiles because they divided the 13 isolates into 3 main 
groups, and differentiated AG18, AG13 and AG14 into 
different vpma subtypes. Strains of pulsotype A had 
vpma profiles similar to that of PG2 except for an add- 
itional hybridizing band in A3F that suggested the pres- 
ence of an additional vpma gene in the Spanish strains 
in this group. All but AG18 also had a size difference in 
the fragment hybridizing to vpmaU compared to PG2. 
Strains in group B, as typed by PFGE, had an identical 
vpma profile that was clearly different from those of the 
other groups and from PG2 and 5632. Finally, the two 
strains isolated from semen (AG13 and 14) were also 
clustered using this approach and were clearly distinct 
from the other groups with a number of vpma differ- 
ences, including the absence of bands hybridizing to 
vpmaXJ, V, Y or Z specific probes, together with a set of 
additional fragments hybridizing to the A3F probe. This 
suggested that isolates from semen possessed a different 
vpma repertoire from strains isolated from milk. These 
data correlated with those obtained with PFGE and 
VNTR. 



Identification of potential new molecular typing tools 

Based on the results obtained above, we assessed whether 
we could further divide the 3 main groups of Spanish iso- 
lates into subtypes. For this purpose, Southern blot ana- 
lyses were first conducted using a set of 7 probes 
designated AGA/BOV that had been identified previously 
as specific for MA or Mycoplasma bovis but that had also 
been found to be subject to inter-strain variation within a 
species (Additional file 2: Table S2). AGA/BOV profiles 
(Table 1) correlated perfectly with the typing clusters 
obtained by PFGE and VNTR analyses, but had no add- 
itional power of discrimination (D = 0.67). 

Previous in silico comparative analysis indicated that 
5632 and PG2 genomes differ by a number of genes, 
most of which fall into two categories: genes that are 
related to mobile genetic elements such as insertion 
sequences (IS) or integrative conjugative elements (ICE), 
and genes that encode enzymes of restriction modifica- 
tion systems [17]. Based on these data, PCR assays were 
then designed that either targeted the presence of a par- 
ticular gene or the genomic location of particular ele- 
ments (see Additional file 3: Table S3 and Additional file 
4: Figure SI). A total of 20 PCR assays were performed 
on our set of isolates, and those with identical results 
were grouped (Table 1). PCR set 3 contained 6 PCR 
assays that yielded identical results with the isolates: no 
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amplification from isolates in the A and C pulsotypes 
and amplification of a fragment of the expected size with 
isolates in the B and D pulsotypes. Since the assays in 
PCR set 3 mainly targeted mobile genetic elements such 
as IS or ICE, these data indicate that isolates from Yecla, 
Malaga and Lorca lack these mobile elements, clearly 
setting the Jumilla strains apart. Southern blot analyses 
further demonstrated the presence of the ICE only in 
the Jumilla isolates. The results obtained with PCR sets 

3 and 6 are consistent with the grouping of the 13 iso- 
lates into 3 clusters, as with PFGE, VNTR and 
hybridization with AGA/BOV probes assays. Isolate 
AG4 could be distinguished from the others with PCR 
set 5 (4 PCR assays). When combined, the PCR sets (or 
at least a combination of 3: set 5 with set 2 or 3 and set 

4 or 6) were as discriminative as PFGE (D = 0.73). 

Discussion 

In the present study, we assessed the use of validated 
and new approaches for the differentiation of caprine 
MA isolates collected from restricted geographical areas 
with endemic CA. Data showed that although the num- 
ber of isolates included in our study was limited, their 
genetic content was surprisingly diverse when compared 
to the previous situation reported for sheep [11,13,15]. 
Thus, typing tools already validated with MA mostly iso- 
lated from ovine, such as VNTR or PFGE assays 
[9,13,24], were shown to be also useful to discriminate 
MA collected in goats from a common geographical ori- 
gin. This paves the way for addressing more complex 
epidemiological questions in the future. Our data also 
indicate that other molecular markers first used in this 
study would be worth evaluating on a larger scale, as 
they might increase the capacity to discriminate isolates 
with very similar profiles in PFGE or VNTR analyses. 
Where isolates cannot be distinguished by VNTR, the 
new assays offer alternative typing tools. Nevertheless, 
our data show that VNTR analyses using just VNTR17 
is a powerful means of quickly discriminating goat MA 
isolates, even when they have been collected from a 
restricted endemic area. 

The situation in Murcia (Spain) was ideal for this 
study because it has been known to have endemic CA 
for several decades [18,19] and the dominance of a na- 
tive goat breed in almost all herds limits the entry of ani- 
mals from other areas of Spain. Although a larger 
number of isolates now need to be typed, data collected 
thus far is already informative because of the congruence 
of all typing techniques. They point to the circulation of 
two clonal groups, one in Yecla (profile A) and one in 
Jumilla (profile B), each with distinctive molecular 
features. The isolates from Yecla resemble the PG2 type 
strain, while isolates from Jumilla are rich in mobile 
genetic elements that are not found in PG2. This 



suggests a different origin of infection with no transfer 
between these two areas, which are only 28 km apart 
(see Additional file 5: Figure S2). The two semen isolates 
from the Lorca area, 150 kms further South, gave almost 
identical profiles that differed from those collected from 
milk (in Lorca and other areas) as indicated by their 
vpma profiles and the presence of sequences coding for 
restriction modification systems (PCR set 4). Whether 
they represent clonal isolates is not known and there is 
certainly a need for detailed studies to determine 
whether isolates from semen form a true subgroup 
within MA. The prospect of associating a particular typ- 
ing pattern with a particular tissue tropism is appealing 
but speculative at this time. To address this issue, more 
isolates collected in this area from semen and from milk 
or cases of mastitis would have to be analyzed. 

In a recent study, a combination of genome-specific 
DNA probes and VNTRs indicated that CA in the 
Western Pyrenees region over the past 30 years has been 
caused by a unique subtype of MA [15]. This indicates 
that eradication programs implemented in the 1990s in 
this area, while efficient in reducing the prevalence of 
the disease, failed to eradicate the strain, which has re- 
emerged over the last few years, causing significant eco- 
nomic losses. The situation encountered in the endemic 
area in Spain differs in that the genetic diversity of the 
circulating strains is high and that several strains appear 
to be present despite the restricted sampling area. This 
situation is likely to be the result of a long period of 
colonization and evolution, but raises the question of 
whether it is a characteristic of MA in goats or whether 
it may also occur in sheep. In goats, clinical signs are 
rare, and sporadic outbreaks are often related to the 
introduction of infected animals from other herds [4]. It 
has been suggested that the genetic diversity of MA 
might be linked to the importation of animals, with MA 
in self-sufficient countries having less genetic hetero- 
geneity [10]. In the endemic area of Murcia, importation 
of goats is not common practice and thus is unlikely to 
be the reason for the genetic diversity of MA. The oc- 
currence of asymptomatic auricular carriers in goat 
herds, with several mycoplasma species cohabiting in 
the ear canal raises the prospect of genetic exchange, 
which may contribute to the genetic variation seen in 
the Spanish isolates. 

Conclusions 

A number of typing tools are now available for the dif- 
ferentiation of goat MA isolates collected in restricted 
geographical areas with endemic CA. The genetic diver- 
sity of MA observed in infected goats in Spain contrasts 
with that seen in sheep with CA in southern France, 
suggesting that assessment of the disease status in 
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endemic areas may require different approaches, de- 
pending on whether it involves goats or sheep. 

Methods 

Isolates and strains 

A total of 13 field isolates were included in this study. 
They came from 3 main breeding areas (Jumilla, Yecla 
and Lorca) located in the goat-producing region of 
Murcia, Spain (see Additional file 5: Figure S2), a region 
known to have endemic CA [18,19]. They were isolated 
mainly from bulk tank milk (BTM) but two isolates were 
also obtained from the semen of 2 asymptomatic bucks 
(Table 1). One isolate from Malaga, a goat breeding area 
located further southwest, was also included. All these 
strains were analysed with strains PG2 and 5632, which 
were used as reference strains and have both been fully 
sequenced [16,17]. 

Culture and DNA extraction 

The isolates were identified as MA by both growth in- 
hibition tests [25], dot immunoblotting on a filtration 
membrane (MF-dot) [26] and a species-specific PCR 
[11]. All strains were cultured in SP4 broth at 37°C [27]. 
Genomic DNA was extracted with chloroform using 
standard procedures [28]. 

Pulsed field gel electrophoresis (PFGE) 

Agarose plugs containing mycoplasma DNA were pre- 
pared as described previously [29]. They were incubated 
in 0.1% Triton X-100 for 2 hours at 4"C before digestion 
with restriction endonucleases overnight at 37°C using 
50 units of either Smal, Mlu\, BstTl and Kpn\ (Promega, 
Charbonnieres-les-Bains, France). Electrophoresis was 
performed using the CHEF-Mapper system (Bio-Rad, 
Marnes-la-Coquette, France) with the standard \ conca- 
tamer ladder (Bio-Rad) as a molecular size marker. 
Agarose gels (1%) were run in 0.5 X TBE at 14°C, at a 
field strength of 6 V/cm, over 24 h with an included 
angle of 120°, and the pulse time ramped linearly from 
10 to 60 seconds. Gels were stained with ethidium 
bromide and visualized using the Gel-Doc 2000 image 
analysis system (Bio-Rad). 

Variable number of tandem repeats (VNTR) 

VNTR analyses were performed as described previously 
[9,15]. Briefly, PCR assays were conducted in a Master- 
cycler ep-Gradient thermocycler (Eppendorf, Le Pecq, 
France) in a total volume of 25 ul containing 0.4 mM of 
each primer, 1 X PCR buffer (with MgS0 4 ; New England 
Biolabs, Ipswich, MA, USA), 200 mM of each deoxynu- 
cleoside triphosphate and 2 units of Taq DNA polymer- 
ase (New England Biolabs). The PCR thermal program 
consisted of 5 min at 94°C, followed by 30 cycles of 
1 min at 94°C, 45 s at 56°C and 1 min at 72°C, and a 



final extension step for 10 min at 72°C. The PCR ampli- 
fication products were analysed by gel electrophoresis 
on 2.5% (w/v) agarose gels and visualised, after staining 
with ethidium bromide using, a UV transilluminator. 

For each VNTR, one PCR product representative of 
each distinct profile was sequenced. The product 
obtained from the type strain PG2 was also sequenced. 
The PCR products were purified using the QIAquick 
PCR Purification kit (Qiagen, Alameda, CA, USA) and 
then sequenced by IF30, CHU Purpan, University of 
Toulouse (France). Tandem repeats in each sequence 
were identified using the on-line tandem repeat finder 
program at http://tandem.bu.edu/trf/trf.html. The data 
obtained were then compared using Needle at http:// 
www.ebi.ac.uk/Tools/emboss/align/ and ClustalW2 at 
http://www.ebi.ac.uk/Tools/clustalw2/ to determine the 
consensus sequence for each VNTR, and the number of 
repeats in each product. Finally, strains were classified 
according to the size of the fragment and number of 
repeats. 

Southern blots 

Several probes were used to examine genetic variability 
in strains of MA. The repertoires of vpma loci in the 
isolates were determined using 6 individual probes based 
on the PG2 sequence, and the general probe A3F 
[21,22]. Profiles were also obtained using a set of 7 AGA 
and BOV probes, previously shown to discriminate MA 
isolates from diverse origins [11]. In addition, the pres- 
ence and number of ICEs in the genome was determined 
using Digoxigenin-labelled probes for ICEA 563 2 (CDS1, 
CDS5 and CDS22) [11,15,17]. Approximately 1 ug of 
genomic DNA was digested with Hindlll, fcoRI (Pro- 
mega Laboratories), Asel or Alwl (New England Biolabs), 
subjected to 1% agarose gel electrophoresis and trans- 
ferred onto nylon membranes (Roche, Indianapolis, IN, 
USA). Specific probes for each of the six vpma genes in 
PG2 (vpmaU, vpmaV, vpmaW, vpmaX, vpmaY and 
vpmaZ) were labelled with Digoxigenin-ll-dUTP 
(Roche) by PCR, as described previously [21]. The 
probes AGA and BOV used in this study were prepared 
by PCR as described previously [11]. 

In all cases, after a denaturation step of 100°C, probes 
and membranes were incubated in Church buffer [30] 
overnight at 55°C. The membranes were then briefly 
washed in 0.2 x SSC/0.1% SDS at room temperature and 
then in 0.2 x SSC/0.1% SDS for 1 h at 55°C. The oligo- 
nucleotide probe A3F, (5'- AA(A/G)TG(T/C)GG(A/T) 
GG(A/T)AC(A/T)A(A/C)(A/T)(A/G)A-3') [15,22] was 
incubated with the membranes overnight at 45°C and 
two washing steps were then conducted in 6 x SSC/0.1% 
SDS for 10 min each at 45°C. In all cases, hybridized 
probes were detected with alkaline phosphatase- 
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conjugated anti-Dig Fab antibodies and CDPstar 
reagents (Roche). 

Identification of new PCR assays as typing tools 

Several PCR assays (see Additional file 3: Table S3) were 
performed to explore differences between the strains on 
the basis of previously observed differences between 
PG2 and 5632 and on in silico comparison of their gen- 
omic se quences [16,17]. Further specific PCR assays 
were conducted. When necessary, primers were designed 
using the on-line program Primer3Plus at http://www. 
bioinformatics . nl / cgi-bin / primer3plus / primer3plus . cgi 
(see Additional file 6: Table S4). All PCRs were per- 
formed as described for the VNTR analysis by incuba- 
tion for 2 min at 94°C, followed by 30 cycles of 30 s at 
94°C/30 s at the appropriate melting temperature (see 
Additional file 3: Table S3) and 30 s at 72°C, and then by 
a final extension step of 5 min at 72°C. PCR products 
were analysed by electrophoresis in 1% agarose gels. 

Additional files 
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